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Abstract. Quantum mechanics is considered to arise from an underlying classical structure ("hidden variable theory", "sub- 
quantum mechanics"), where quantum fluctuations follow from a physical noise mechanism. The stability of the hydrogen 
ground state can then arise from a balance between Lorentz damping and energy absorption from the noise. Since the damping 
is weak, the ground state phase space density should predominantly be a function of the conserved quantities, energy and 
angular momentum. 

A candidate for this phase space density is constructed for ground state of the relativistic hydrogen problem of a spinless 
particle. The first excited states and their spherical harmonics are also considered in this framework. 

The analytic expression of the ground state energy can be reproduced, provided averages of certain products are replaced 
by products of averages. This analysis puts forward that quantum mechanics may arise from an underlying classical level as a 
slow variable theory, where each new quantum operator relates to a new, well separated time interval. 

INTRODUCTION 

In this centenary of Einstein's Annus Mirabilis it is appropriate to reconsider the foundations of physical theories. 
We shall, in particular, be interested in a possible bridge between the three fields that Einstein pushed so much in 
1905: brownian motion, special relativity and quantum mechanics. Our focus is to investigate the possibility that 
quantum mechanics arises as a statistical theory from a certain special relativistic theory, with specific stochastic 
forces constituting the notorious "quantum noise". Conceptually this will be a "unified field theory", which Einstein 
tried to formulate in the second part of his life. Without fully specifying this theory, we shall generalize squares of 
quantum wavefunctions to phase space densities and see what it leads to. 

In order to motivate our investigations, let us notice that the status of quantum mechanics is still judged in two 
basically opposite ways. The empirical view stresses that quantum mechanics, with its many postulates such as the 
superposition principle, the uncertainty principle, the exclusion principle, the Born rule, the collapse in measurements, 
is our best theory of nature, proven to work in all cases where it was tested, ranging from solid state physics 
and quantum chemistry to high energy physics and the early Universe. Most physicists even believe that quantum 
mechanics is the ultimate theory of nature, that must be unified with gravity, one popular route being string theory, 
which adopts the additional postulates of supersymmetry and extra dimensions. 

On a parallel track, there has always remained the fundamental question: what is quantum mechanics actually 
standing for? Does nature have a reality: do particles exist as physical entities, possessing a position and a momentum; 
is the moon there when nobody looks? What is the physical mechanism behind the notorious "quantum fluctuations" 
and uncertainty relations? What causes interference effects and the particle-wave duality? Does each particle have its 
own wavefunction (standard Copenhagen interpretation), what does a "collapse of the wavefuction" mean physically 
' ? For decades these questions have suffered from the lack of progress, in particular in the quantum measurement 
problem, the only undisputed point of contact between quantum mechanics and nature. Very different interpretations 
have emerged [2]: the Everett- Wheeler relative state or multi-universe picture (collapse does not occur, but rather a 
branching); the mind-body problem (the registration of the measurement in our mind or a machine matters) (Wigner, 
see [2]), the wavefunction codes our state of knowledge [3] or our state of belief [4]. 



Nothing measurable according to Zbinden et al [1]. 
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At present a variety of different approaches support the possibihty that at the quantum level, nature has a 'natural' 
(i.e. more classical than long supposed) behavior. We may mention: arguments against the conclusion that nature 
should be non-local because of violations of Bell inequahties [5, 6, 7]; the demonstration that quantum probabilities 
including interference effects are compatible with classical probability theory [8]; the long effort of Stochastic 
Electrodynamics [9] and its later version called Linear Stochastic Electrodynamics [9, 10]; the demonstration of 
entanglement in classical brownian motion [11]; the description of an electron with its spin as a soliton (over-extremal 
Kerr-Newman black hole) in electro-gravity [12, 13]; the many arguments in favor of the statistical interpretation 
[14, 15, 16, 17, 18, 19]; explicit solutions for the quantum measurement problem [18, 19]. 

Let us go in some detail about the measurement process. By taking into account an explicit model for the apparatus, 
the quantum measurement problem was recently solved in our group in two situations: in the first the position of a 
boson is measured by an apparatus that consists of many bosons. The apparatus starts close to an (ideal) Bose-Einstein 
phase transition and is driven into it by the measurement [18]. In the second model the z-component of a spin-^ is 
measured by an apparatus consisting of many spins- j [19] that start out as a paramagnet and are driven into the up 
or down ferromagnetic phase, according to the sign of the measured spin component. In both cases the apparatus also 
consists of a heat bath, which originates from other degrees of freedom, in the magnetic case the lattice vibrations. 

Quantum mechanics being a probabilistic theory, it should come as no surprise that what could be calculated was the 
statistical distribution of outcomes for an ensemble of measurements on an ensemble of identically prepared systems or 
identical preparations of a single system. Technically this comes in because the whole exercise is to determine the post- 
measurement density matrix given its pre-measurement form. The process goes in two steps: when coupling the system 
to the apparatus, there appears a very fast collapse of the off-diagonal terms (the death of Schrodinger cats). This 
"collapse of the wavepacket" is a unitary evolution in the large Hilbert space of system and pointer degrees of freedom 
of the apparatus (thus not involving the bath), very similar to the process of dephasing in NMR physics. The collapse 
time depends as an inverse power on the number of particles in the apparatus, so it is strongly microscopic when 
the apparatus is macroscopic. On a later stage the effect of the bath is to make this collapse definite by decoherence. 
Schrodinger cat terms thus disappear automatically, because of the macroscopic size of the apparatus. The second 
step of the measurement is its registration, which involves the diagonal elements of the density matrix and appears 
to have classical features [20]. Indeed, for the magnetic model, the solution of this part of the problem is just what 
one might have anticipated from purely classical thermodynamic reasoning [20]. In particular, a bath is needed simply 
because energy has to be dumped when the macroscopic pointer variable (the magnetization of the apparatus) goes 
from its initial metastable paramagnetic value to its up or down ferromagnetic value. Also in cloud chambers, the 
energy released upon creating bubbles in an oversaturate liquid has to be transfered to the envirormient and this very 
mechanism allows the selection of circumstances where these bubbles have a desired average size. 

From the recent solutions of two realistic models for the quantum measurement process, it became clear that the 
statistical interpretation of quantum mechanics allows a simple interpretation of this process during and after the 
measurement. To mention one problem that is immediately solved in this way: the fundamental question how quantum 
mechanics can describe collapse of the measured operator has the answer: quantum mechanics describes only the 
statistics of those collapses; it is not the proper theory to describe an individual collapse. Spontaneous collapse models 
(GRWP theories), that try to treat individual collapses from a slightly extended quantum mechanics, attempt to heal a 
problem that has its roots at a more fundamental level. Indeed, for a two slit experiment quantum mechanics does not 
describe any individual spot on the photographic plate, but only the density of many spots. 

Thus, we conclude that Quantum mechanics is not a theory of nature itself, but a theory for the statistics of outcomes 
of experiments. This situation is very similar to the one of thermodynamics, where the underlying level of, say a 
Lennard-Jones gas, is Newton's equations of motion. The main difference is, however, that experiments have to be 
made with apparati which themselves are always macroscopic and therefore influence the microscopic dynamics. 
Since it is the task of physicists to describe nature, the search for a sub-quantum theory (hidden-variables theory) is, 
in our eyes, completely justified. The program to search this may contain the following steps 

1) For point particles, show compatibility between classical phase space density structures and quantum structures. 

2) Find the underlying stochastic theory and isolate the stochastic forces which cause the "quantum fluctuations". 
Show compatibiUty with quantum mechanics. 

3) Find soliton-like solutions for elementary particles such as the electron. An interesting attempt sees an electron 
as a Kerr-Newmann black hole, that is, as a charged ring with diameter given by the Compton wavelength [12, 13]. 
On the ring a current runs with the speed of light, causing a spin, with ^ = 2, as in the Dirac theory for the electron. In 
these theories the parameters must be fixed to the physical values by adding terms to the Einstein-Maxwell theory. A 
possibility is to couple torsion with the electromagnetic field tensor [21], which might fix charge and/or spin, though 
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not yet the mass. Within this framework the PauU principle might arise as an energetic constraint on the behavior in 
the far field region, and Einstein photons as massless spin-one solitons. Let us notice that indistinguishably is a natural 
aspect of sohtons that wind up in a chaotic motion. 

4) Show that the consistency still works for the "stochastic soliton mechanics" of sohtons in the random electro- 
magnetic background, with "photons" appearing as classical EM waves and as "Einstein" solitons. 

Clearly, the program may run into contradictions at many points: apart from internal consistency, compatibility with 
a huge amount of experiments is required. But if Einstein was right after all, then there should be some road to the 
underlying structure, and this one seems to us the simplest, so it deserves to be investigated further. 

In this work we shall restrict ourselves to point 1) of the program and observe that already this one is more 
complicated than sketched here. 



MOTIVATION FROM STOCHASTIC ELECTRODYNAMICS 

Quantum mechanics has both a particle and wave nature. In a subquantum approach, it is clear how to think of 
the particles, to start they are just considered as points, so the next question: what are the waves made of. In the 
simplest approach, one assumes that the waves are due to fields that we know already. The first candidate then is 
Electromagnetism. One has then arrived at the so-called Stochastic Electrodynamics (SED), where the zero-point 
fluctuations known as a slogan in quantum mechanics are attributed to a physical entity, namely classical random EM 
fields [9]. To get the proper cormection, one then has to choose the energy content of each mode as jh(0, in which the 
prefactor h is just a system parameter, chosen to coincide with Planck's constant. The related spectral density is just 
the zero point part of the Planck spectrum, p (co) = hco^ /{2kV). This ©^-dependence is precisely the one consistent 
with Lorentz invariance and with thermodynamics. 

The theory of SED has a long history of both success and failures. SED has explained the harmonic oscillator and, 
with it, electric dipole harmonic oscillators, the Casimir and Casimir/Polder effects, and the Unruh effect. 

Away from harmonic systems, SED has a long record of failures. First, a connection with the Schrodinger equation 
has never been found. Second, the role for excited states has been discussed only in the margin. It should be admitted 
that, whereas Schrodinger dynamics at first neglects the decay of excited states, the structure of SED is that of time- 
dependent perturbation theory, where the full time-dependence is taken into account from the start. Third, a long list 
of non-hnear problems has been studied, where results different from quantum mechanics emerge, see e.g. [22]. 

The most obvious example wiU be the subject of this study: the ground state of the hydrogen problem. This is a 
rather clean problem that should be understood before conclusions are drawn about vahdity or invahdity and about 
multi-level atoms. Indeed, in the hydrogen problem no internal structure should be expected that complicates the 
orbiting of the point-electron in the central force field of the infinitely heavy nucleus. Theoretical results so far are 
meager, predicting stabihty radius of the orbit equal to infinity [9]. Surprisingly, a numerical analysis by Zhou and 
Cole supports an SED explanation of the quantum ground state density [23]. 

That the hydrogen problem should have a stable solution within SED, was put forward long ago: when the electron is 
far from the nucleus, its orbital frequency is low and so should be the relevant frequencies of the stochastic force. Due 
to the or' spectrum, there should essentially be no random force at all. What remains is the Lorentz damping that brings 
the electron closer to the nucleus. For small distances, on the other hand, the relevant frequencies should be large, and 
the spectrum predicts a strong kicking due to the random forces, preventing electron to fall into the nucleus due to 
Lorentz damping alone. There thus appears a stable atomic state. Is this the one we know from quantum mechanics? 
So far attempts to show this have failed, including those by ourselves [21]. For this reason it would be interesting to 
corroborate the numerical analysis of [23]. 



QUANTUM MECHANICS OF THE RELATIVISTIC HYDROGEN ATOM 

The Hamiltonian for a spinless boson with rest mass m and charge e, in the presence a central charge —Ze, reads 

// = -y/mV-hpV = Ymc^ (1.1) 

r r 
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We express r, p and time in atomic units, 



^^r, 1^., P^^p (1.2) 
ao To mao 



with Bohr radius and Bohr time 



^ ^ ^ ^5.2910-"., ro ^ = ^2.418 10"". (1.3) 

with numerical values for the electron. This scaling also implies 

E L 

— ^E, -^L, (1.4) 
mc^ n 



and brings the scaled form 

H = ^/l + ^^^-^^=Y-^^ (1.5) 



r r 

Let us denote quantum operators by a subscript op. The Schrodinger equation reads 



{Hop - E)\xif) = ( a/i + «'^Pop -^-E)\xif)=0 (1.6) 



this is equivalent to 

^2-7-2 77 1 172 



,1 2 a^Z-" £ M , , 

(2P^p-^-;-)lv^) = -^lv^) (1-7) 



In the position basis this reads 



When seating r ^ r/E,p ^ Ep, and using thatL^^ has eigenvalues /(/ + 1) with non-negative integer /, this becomes 

,1^2 1 l(l + i)-a^z^ r 1 

+ J ' _y,= _ (1.9) 

2 r 2r^ r /? 



so —1/R is the associated non-relativistic energy. It is customary define the effective angular momentum value £/ by 

eiiei + l) = l{l + l)-a^z\ ii = -i + y{2l + l)^-4a^Z^^l-^ (1.10) 
Now the whole problem can be solved by analogy with the non-relativistic situation. The ground state wavefunction is 

The so called Yrast states have (n,/) = {n,n — 1) and wavefunctions 

= C/'e-''^'"/(^'+i) (1.12) 

They have energy 

En,i = {1 + f^' « 1 - ^ - «'2^[375|tTT - A ] (1-13) 

a formula which is valid for all levels, that is to say, forn = 1 , 2, • • and < / < n — 1 . 

Let us mention that when the spin of the election is taken into account, one has to replace / — > j = l±j. 
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STATISTICAL MECHANICS OF THE RELATIVISTIC COULOMB PROBLEM 



The Schrodinger equation works so well for the hydrogen atom, because corrections arising from the coupling to the 
electromagnetic field are weak. They lead to the Lamb shift, which is of order a^Z"* In aZ in the present units. So it is 
weaker than the fine structure, represented by the af'Z'^ term in (1.13). The very same mechanism makes the lifetimes 
of excited states much larger than the Bohr time. 

To focus our Une of reasoning, let us consider the theory of stochastic electrodynamics (SED). Here random 
electromagnetic fields are assumed to exist. In this theory a Fokker-Planck approach may be formulated for the 
statistics of the random orbits (r(/),p(f)) [9, 21]. The phase space density will satisfy an equation of the form 

/'(r,p,0 = (5(r(0-r)5(p(/)-p)) (2.1) 

will satisfy 

dtP= —££P + Lorentz damping + diffusion terms (2.2) 
with drift given by the LiouvilUan 

^=ipVr-4'^P (2.3) 

Since damping and diffusion terms are of order w" , the stationary distribution should, to this order, be a function of the 
conserved quantities, angular momentum L and energy E. This remains vaUd if one goes beyond the Fokker-Planck 
equation and for other stochastic theories. So the weakness of the Lamb shift finally allows us to consider equihbrium 
structures of a class of stochastic theories, without going to specific details. 
The ground state should follow from the dynamical problem H = j— a}Z^ / r = E or, equivalently, 

1 2 «^Z^ E _ E^-\ 

2^ ~2i^ 7~2^ ^ ' 

There wiU also be effects of damping and noise, but they enter here only at order a^, and should explain the Lamb 
shift. We shall neglect those effects in the present paper. If we define R by 



the dynanaical problem becomes 



1 2 «2z2 E E 



2" 2r^ r R 

Close to r = and for pr = Oit is seen that pj_ = aZ/ r, implying that L > aZ and R>r. Let us denote 



0) = y 1 - (o^L^ =L^- a^Z^ (2.7) 

We shall take spherical coordinates in r space. In the frame along r = r(0,0, 1) we introduce for p the angles and 
V with 0<n<n and < V < 27r, 

P — {P-LiTPi-2,Pr) — /?(sinjUcosV,sinjUsinv,cos/i). (2.8) 

In special relativity the invariant volume elements in coordinate and momentum space are 

dVr = d^r = drded^r^ sin 0, dVp = ^ (2.9) 

We have, using that LdL = C0Ld{(0L), 

dprdvdp±p± dprdvLdL dprdv (0Ld{(0L) 
uVp = = 7^ = ~ yZ.ly)) 
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We go to the new variables n and R by the transfomation 



and have 



(0L = r\j2E{^-^) Sinn, pr = J 2E{-^-j) cos ^ (2.11) 



, , ,„ ^ ,E rdudRvE ^ a^Z^ , , 

dprdi(0L) = rdudR^( ) = , Y=E-\ (2.12) 



This implies the volume element in momentum space 



dV„ = , , \ smu ^(E) = ^ ^ - ' — (2.13) 

2R^<$>{E) \ Er ER ^ ^ ' 2E^ 

Notice that the r-dependence through 7 has dropped out and that ^ = j non-relativistically. If the phase space density 
has the form ^{E,L) = (oLR^(t>{E) exp{—aR), then we have 

dVp » = d\ldvdR {R - r)e''''^sm^ (2.14) 



The momentum integrals, that is to say, the integrals over fj,, v and R>r, then yield 

n2 



^ e-"' (2.15) 



which for a = 2 is just the square of the non-relativistic ground state wavefunction. Furthermore, since 

()p-l}R 



2E 



= r{R-r)sm^H, (2.16) 



powers of r in the squared wave functions can be generated by powers of this ratio in the phase space density. So we 
may choose the following combination of the conserved quantities E, R{E) and L for ground state (1,0) and, more 
generally, the Yrast states (n,/) = (n,n — 1) with n = 1,2, • • • 

^{R,L) = C(OLR^<^{E) (2.17) 
Fixing C this leads to our final proposition 

rlu 96+8*', p5+4C, 

This form is very appealing, though perhaps not unique. The angular integrals yield 
and the fuU momentum average brings the proper squares of the wavefunctions, 

yp''^''^"Vv " (i+^/r+2^'r(2+2^,) ' ™ 



On the other hand, the spatial average yields 



5+4f, 



dVr, Hdrr^^^dR^^ , "" ^ sin2+«' ^ (2.21) 

^io n 2n (i+^,)6+4£,r(| + 2^/)r(|+2^,) 

Classical Phase Space Density for the Relativistic Hydrogen Atom February 1 , 2008 6 



The angular integral brings 

/ an — AH 

(l+£,)f'+4</r(4 + 4£;) 



fdr? [ 
Jo Jn,v 



23+41'/ £■5+4 



dVp^ = dR ^^4^^^^^±^L___/?4+4^,g-2i;£„,/(i+^,) (2.22) 



This impUes the expectation values 

1 ^ Eni ^ ,2 j_ 2E^ (—) = ^ (2 23) 

(l+^;)2 ^R'' (1+^,)3(1+2^,)' ^/?2/ (1+^^)3(3+4^^) ^ • ^ 

Eq. (2.5) then implies, to order a^, the classical expectation value 

. = 1 - a^Z^(^) + -^a^z\-,) = 1 - ^ - V,4(s,2_J+i) = + ^^^.^^ (2.24) 
This is not the quantum mechanical result, so the approach fails at order a^, that is to say, for the fine structure. 



NEGLECTING CORRELATIONS 

It is easily verified that the proper a'^ term of the quantum eigenenergies (1.13) emerges is we replace the term {R^'^) 
of Eq. (2.24) by {R~^)^, thus neglecting the correlator {R~^) — {R~^)^- This approach appears to work to infinite order. 



Yrast states 

Surprisingly, the quantum mechanical result does emerge if, after expanding the square root (2.5) in powers of a, 
we replace in the successive terms the averages of products by products of averages. The result can then be resunnmed. 



From the expression (1.13) one may derive the relation 

a^Z^=^—^{n-l + i,f (3.2) 
For the Yrast states we then observe agreement with the quantum eigenvalues for the energy. 



to all orders in a. Physically such infinite precision of fortuitous, since there are Lamb shift corrections of order a^. 

By working with expression (2.4), we have another approach, which acts as a consistency test. If we separate the 
averages and calculate the average of E as above, we find 

= ((2£(1 - 1) cos2;u)) = 2{{E)) ((1 - 1) cosV) = (3.4) 

and 

,0)21.2 E2 (3+4^0 
^2 ' 

These results bring 



^ -2 / /Id J- ^ ^ 
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Combining both results brings simply 



L2 , 1 



((P')) = ((P?)) + ((^)) = ^YT^ (3.7) 

When averaging Eq. (2.4), we again expand in powers of a}"!? and replace repeatedly averages of products by products 
of averages. We end up with the alternative expression 



mUci = sjl + a^^m)- ((^)) = ^1 + ^ - ^ (3.8) 

The leading correction of this square root is stronger than the one of order a'^Z'^ of the expression (3.1). But by using 
eq. (3.2) is easy to see that 

mUci = \l^-^-^=Eni (3.9) 
confirming consistency of the approach, provided it includes the replacement {{r~'^)) {f'^)^ done in Eq. (3.6). 



The 2s state 

The Yrast states could be special, because their spatial wavefunctions consist of a single term. We now investigate 
a non- Yrast state, the 2s state, which has and additional polynomial in its wavefunction: 

= { (2+4)(4+34?3/2+4r(2+24) } (i+4)V4+34)^''K"V?+1?^) ^^'^^^ 

The function /?2o(^) is orthogonal to ^io(^)> as it should in view of the scaled Schrodinger equation (1.9). The 
probability density P2o{>') = ^2o('') derive from a phase space density, 

^-"^^'"'^^^^^^ (2 + 4).^r(6 + 44)(4 + 34)^/^^^^o -^'^(^)[-^) -K-7?mJ 

2 2 ^ 2 2 \ ^ 

(1 +4)2(4 + 34)2(3 +44) (5 +44) - 8(1 + 4) (4 + 34) (5 + 4^) + 16 [^j^) (3.11) 
which after integrating over the angles n and V reads 

X [(1 + iof (3 + 24) (4 + 34)^ - 2(1 + 4) (3 + 24) (4 + 34)K^ - '") + ^r\R - rf] 

and next 

X [(1 + 4)^(4 + 34)^ (5 + 44) (7 + 4 4) - 2(1 + 4) (4 + 34) (7 + 44)/?^ +/?^ ] 

This brings the averages 

^7^" (2 + 4)^4 + 34' 2(2 + 4) v?+34' ^^'^"^^ 
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which is just what is needed, because we get using a^Z^ = — (^q + 1 ), 



1, 



{{E))20;cl = \ll + a^Z^{-r-{ ^ 



) = 



: + 



2 + i 



We can further verify that 



((p2))^2((£))[(i)-(l)]+a2z2(i)2 



1 



(2 + 4)2 

implying that also the alternative calculation leads to the quantum mechanics result, 

4(4 + 1) 



= E2o (3.15) 



(3.16) 



+ 



2 + in 



V4 + 34 

2 + 4 ' (2 + 4)\/4 + 34~ \/4 + 34 



It appears that the approach is non-unique. The density /?2o('') obtained from 



'^0 



26+44(1+^,,) 



(2 + 4)r2(3 + Ito) (4 + 34)^/2+2^0 



2(0 



2E 



exp 



2R 



( 1 + ^o) (4 + 34) (7 + 84) -8(1+4) V4 + 34/? + 2R^ 



2nB{^+2£o,^-) 



V4 + 34 
1 -1 



(3.17) 



(3.18) 



This form has the same expectation for (1/ r) and (l/R), and thus for energy. The two forms can thus be mixed, which 
leads to a one-parameter family of solutions; the proper choice cannot be fixed in the present approach. By viewing 
such mixings as a way to fix the term of the square of the polynomial of R20, one might wonder whether the r term 
can be changed with a second independent parameter. It has been checked, however, that that would modify ( 1 /R) . 



PHASE SPACE FORMS FOR SQUARES OF SPHERICAL HARMONICS 

Previous expressions R^q for s = states are still multiplied by a factor Fqq =1/ (47r), which drops out after the 
integration over the angles 6 and <j) of the coordinate vector. Let us denote the generalization of Yf^Yii^i by . 
Then we can likewise multiply our results for and ^20 by ^^q" = 1 /(4;r) and cancel it by the angular integrals. 

For the ^21 distribution we only considered the radial part; there is still the question how to deal with the factor Y^^ 
and the same question arises for the higher / states at n > 2. According to the quantum solution, such diagonal terms 
should have no time-dependence, so also Fj^^ is be generalized to a function depending on conserved quantities only. 

Let us notice that in the frame along r we can introduce angles fj. and V with < n < n and < v < 2n, 

r = r(0,0,l), p =/7(sin/i cosV,sin/i sinV,cos/l) = (/7_Li,/7_L2,/>r), L =/7rsin/i (— sinv,cosv,0) (4.1) 
In the laboratory frame this impUes 

r = r(sin 0cos0,sin0sin0,cos0), L = prsin/i L, 

L = (— cos cos ^ sin V — sin ^ cos V, — cos sin ^ sin V + cos <j) cos V , sin sin v) (4.2) 
Specifying the z-axis as our preferred axis, we have = sin sin v and we define 

In general with n a fixed unit vector this would read 

= A [(n . t)2 - „ . t] , ^jio = A [1 - 2(n . t)2] , ^,11 = A [(n . t)2 + „ . t] , (4.4) 

The odd terms are chosen linear, as this seems most natural to us. To settie their form uniquely is, however, not possible 
in the present approach, since for that aim the underlying physics should be specified further. 
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The p-integration produces an integral over v, which brings the desired results 
Moreover, it holds that 

= ^sin20, (m = -l,0,l) (4.6) 

If we integrate this over all angles 6, (j) and make the identification L^gp <-> we reproduce the quantum result 

Lz^im-im= J smdddd^ j =m= {Y,J^L,op\Yx,n), {m = -1,0, 1) (4.7) 

and also the related results 



= = {Ylm\Uop\Yl,n), = = {Ylm\Ly „p\Yx,n) ■ (4.8) 

We conclude that at least the lowest non-trivial square of spherical harmonics can be properly generahzed. 



COMPARISON WITH THE WIGNER FUNCTION 



The Wigner function of the ground state wavefunction reads in the non-relativistic Umit [24, 25] 

2^n /■ 3 

{iKfJ ^[l + (p+lq)2]2[l + (p_lq)2]2 



As opposed to our phase space density, this function is partly negative, a fact commonly considered as a key argument 
against any underlying quasi-classical reality. Its marginal reads 

The non-relativistic limit of our phase space density (2.17) reads for the ground state 

^(r,p) = ^LR^e-^^ (5.3) 
Taking spherical coordinates along p yields L = sin 9 and using — 1/r = — 1/^ one derives for its marginal 

^(p) = ^ [^'"'drr^Rh-^'' = ^ HdR f^e'^* (5.4) 

71 Jo n Jo il + hp^R)^ 



(l + ip2/?)5 

Whereas (5.2) is finite at /? = and decays as p^^ for large p, Eq. (5.4) is linear at /? = and decays as for large p. 
Figure 5.1 compares these results. It is seen that the difference is mainly at small p, where because of the phase space 
volume little weight is located. The tails are already small when they start to deviate from each other. 



DISCUSSION 

In this paper we have considered the question whether the quantum probability densities (the squares of quantum 
wavefunctions) of the relativistic hydrogen atom, in the absence of electron spin, may arise from the phase space 
distribution in an underlying classical stochastic theory ("sub-quantum theory" or "hidden variable theory"). Because 
of the small value of the fine structure constant, this is a weak damping problem, for which it is generally known 
that the stationary distribution is determined by conserved quantities (energy and angular momentum) alone. Without 
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FIGURE 5.1. Reduced phase space density ^{p) of the non-relativistic hydrogen ground state, Eq. (5.4), (lower curve), and 
reduced Wigner density W(p), Eq. (5.2), (upper curve), as function of the momentum p, in atomic units. They are both normalized 
to unity. 



specifying the underlying stochastic theory, we have shown by construction that such a scheme can work for the 
considered cases. The relativistic version of the hydrogen problem has been chosen in order to have the fine structure 
constant as a parameter, and the generalization works to all orders (be it that from the outset, all Lamb shift effects are 
neglected, in our approach as well as in Schrodinger dynamics). 

Our resulting phase space density ^(r,p) does not coincide with the Wigner function W(r,p). For the hydrogen 
ground state the latter can be expressed as derivatives of an integral [24, 25]. In contrast, for our ground state ^^(r,p) 
is explicit, it contains powers and an exponential. Although we did not prove its uniquness, it looks simple and has been 
generalized for the whole family of Yrast states. It is well known that the Wigner function for the hydrogen ground 
state is partly negative, and this has long been considered as a no-go argument against any classical-type subquantum 
structure. However, our ^(r,p) is non-negative for the ground state. Does this solve that paradox? Yes, because our 
momentum stands for the instantaneous momentum of the particle. The momentum of the Wigner function is defined 
in via off-diagonal elements of the quantum density matrix, and should be considered as a coarse grained momentum, 
just as would happen in brownian motion [26, 1 1]. For such a coarse grained momentum positivity of the density is 
not obvious, probably not needed, and certainly absent. 

From a physical point of view, there are many reasons to accept that particles are to good approximation point 
particles: traces in cloud chambers, images of non-flat crystal surfaces obtained by scanning tunnel microscopy, 
trapping of ions in magnetic fields over periods exceeding a year. If, one the subquantum level they just have a position 
and a speed, the phase space density cannot be negative. The problem of negative probabilities is solved for our ground 
state, but still a potential danger for excited states. For the 2p state the generalized spherical harmonics are partly 
negative, while for the 2* state the radial phase space density is partly negative. But in the theory of classical brownian 
motion the excited states are also partly negative, while the total probability density will always be non-negative. 

Within Stochastic Electrodynamics the 2s state can be populated by a two-photon excitation from an ensemble of 
ground state atoms, but it cannot happen in such a manner that the total phase space probability is negative anywhere. 
For instance, one could imagine a large population of the Is state and a partial one of the 2s state, with not- too-large 
values of L and R in the 2s state, such that the total probability is still positive everywhere. Afterall, at the classical level 
each realization consists of an electron in some orbit around the nucleus in the presence of fluctuations, and a "state" 
merely describes an ensemble of such systems. The 2p phase space densities at m fixed are partly negative, so they 
cannot exist on their own. However, the 'transversally' polarized mixture ^2i{E,L){^^^Si +^n) is non-negative, 
while the unpolarized mixture ^2i{E ,L)(^^^S^ +^o'' + ^ii') = 3 ^21 {E,L)/{4n) is strictly positive. They can in 
principle exist without ground state buffer. The odd parts may differ from the assumed linear ones of (4.4). When they 
are cubic, bofli ^iL"/ and can be positive, but not i^iQ- The quantum operator L^op may then be identified with 
the classical times a constant. 

Given the phase space distribution, we have considered the average energy. We have observed that the quantum 
energy eigenvalues are reproduced provided certain correlations are neglected, after which averages of products reduce 
to products of averages. Physically this occurs when these averages are taken at very different times, such that their 
fluctuations are uncorrelated. This hints at the possibifity that the operator structure of quantum theory arises from the 
physics of timescale separated observables in an underlying stochastic theory. The property that repeatedly applying a 
quantum operator on one of its eigenstates produces just powers of the related eigenvalue has a natural interpretation 
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in terms of temporal separation of averages, as is long understood in the theory Stochastic Electrodynamics [9]. It is 
remarkable that such a separation appears necessary here to estabUsh the proper mapping of energies between the class 
of stochastic theories and quantum mechanics. 

In atomic physics spectral emissions are due to transition from an excited state to a lower one. The 2p-state has 
a lifetime given by oc^^Tq muhiplied by phase space factors, the final result is of order of magnitude of the inverse 
Lamb shift frequency ~ Ins. So about a^^ Keppler orbits are traversed before the emission takes place. This large 
number of revolutions, together with the gradual influence of the noise that slowly changes the orbits, may explain the 
sharpness of spectral lines within our picture. If only a few orbits are traversed, the sharpness of the spectral line will 
not be achieved, since the spread in energies of Keppler orbits is large. But this just provides a physical picture for the 
energy-time uncertainty relation know from quantum mechanics. 

There are many remaining and often worrying issues. An essential item is the role of dynamics, such as the ls-2s 
and \s-2p interference terms. Likewise, there is the task to include spin. It is good to realize that our present radial 
distributions are functions of L^, which remains conserved in the presence of spin. 

For scattering and ionization experiments on of e.g. electrons on e.g. atomic hydrogen [27] it is known that the 
energy loss of the incoming electron can only take one of the discrete values needed to bring the ground state electron 
to the one of the excited levels (Lyman lines). In the recent work [28] transitions to n = 2, 3 and 4 are clearly observed, 
while the higher states form a quasi continuum, that matches the ionization continuum. Whether these discrete features 
can be described in the present formahsm is a question related to the existence of excited states, and the timescale on 
which they express themselves. What we have discussed is a number of features that should be respected in our class of 
the underlying theories. Whether the proposed type of ground state and excited states exists and what happens during 
electron scattering events, are questions that cannot be answered before specializing the precise type of stochastic 
noise. 

Like any new theory, the present approach poses more questions than it solves. For any sub-quantum theory a vast 
amount of experimental constraints is to be respected. But we are convinced that no individual quantum experiment 
is described by any current theory. So if Einstein was right after all, quantum mechanics is incomplete and then there 
should be some path to an underlying, more complete theory. We hope that the present work plays a role in uncovering 
it. 
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